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■ Aims. We used data from the 58 month long, continuous Swi/f/Burst Alert Telescope (BAT) observations of the five brightest Seyfert 
galaxies at hard X~rays, to study their flux and spectral variability in the 20-100 keV energy band. The column density in these 
objects is less than 10^"* cm"^, which implies that the SwiftfBKI data allow us to study the "true" variability of the central source. 

, Methods. We used 2-day binned light curves in the (20-50) and (50-100) keV bands to estimate their fractional variability amplitude, 

and the same band 20-day binned light curves to compute hardness ratios and construct "colour-flux" diagrams. We also considered a 
thermal Comptonization model, together with a reflection component with constant flux, and produced model "colour-flux" diagrams, 
(— I ' assuming realistic variations of the model parameter values, which we then compared with the observed diagrams. 

Results. All objects show significant variations, with an amplitude which is similar to the AGN variability amplitude at energies 
I . below 10 keV. We found evidence for an anti-correlation between variability amplitude and black hole mass. The light curves in both 

O ' bands are well correlated, with no significant delays on time scales as short as 2 days. NGC 4151 and NGC 21 10 do not show spectral 

variability, but we found a significant anti-correlation between hardness ratios and source flux in NGC 4388 (and NGC 4945, IC 4329, 
C/3 , to a lesser extent). This "softer when brighter" behaviour is similar to what has been observed at energies below 10 keV, and cannot 

■ be explained if the continuum varies only in flux; the intrinsic shape should also steepen with increasing flux. 
Conclusions. The presence of significant flux variations indicate that the central source in these objects is intrinsically variable on 
time scales as short as ~ 1-2 days. The intrinsic slope of the continuum varies with the flux (at least in NGC 4388). The positive 
"spectral slope-flux" correlation can be explained if the temperature of the hot corona decreases with increasing flux. The lack of 
spectral variations in two objects, could be due to the fact that they may operate in a different "state", as their accretion rate is less 
than 1% of the Eddington limit (significantly smaller than the rate of the other three objects in the sample). 
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1. Introduction the corona can also illuminate the disc, producing the Fe Kq, 

line (6.4-6.97 keV) and t he "reflection hump" at higher energies 
T-H , Active Galactic Nuclei (AGN) are the most prominent and per- ^George & Fabian. 199 Ih . Reflection of the central radiation on 
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sistent X-ray sources in the sky. The current paradigm for the molecular torus often le ads to an additional, p ossibly domi- 
central source m these objects postulates the presence of a cen- ^^^^^ reflection component (iGhisellini etal.. 19941) . 
tral black hole (BH) with a mass of 10 - 10 Mq, and a geo- 



5— I metrically thin, optically thick accretion disc that may extend to 
5^ the innermost stable circular orbit around the black hole. Current The AGN X-ray emission is strongly variable, even on time 
unification schemes for radio-quiet AGN also postulate the pres- scales as short as a few hundred seconds, both in flux and in 
ence of an obscuring molecular torus around the central source the shape of the observed energy spectrum. In most cases, radio- 
which prevents a direct view of the continuum and broad-line quiet objects show a positive spectral slope-flux relation, i.e. a 
region in Seyfert 2 galaxies. The accretion disc is thought to softer-when-brighter behaviour, in the 2-10 keV band bofli on 
be responsible for the broad, quasi-thermal emission component short and long time scales (see e.g. Sobolewska & Papadakis 
in the optical-UV spectrum of AGN (the so-called "Big Blue 2009 and references therein). From a theoretical perspective. 
Bump"). At energies above ~2 keV, a power-law like compo- intrinsic slope variations are expected (Haardt, Maraschi & 
nent is observed in radio-quiet AGN. This is atti'ibuted to emis- Ghiseflini 1997; Coppi 1999; Beloborodov 1999). For exam- 
sion by a hot corona (T ~ 10^ - Iff K) overlying the thin pie, in the case of thermal Comptonization models, valuations 
disc. The corona up-Comptonizes the disc soft photons to pro- in the UV/soft X-ray photons can alfect the slope of the X-ray 
duce the hard (E~ 2- 200 keV) X-ray emission. X-rays from specti'um and result in "softer" X-ray spectra when the source is 
bright. In fact such variations may have already been observed 
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On the other hand, it has been argued that this apparent spec- 
tral steepening with increasing source flux in the 2-10 keV band 
does not necessarily represent an intrinsic change in the spectral 
slope. These variations can be explained if the power-law contin- 
uum is variable in flux, but not in shape, and the reflection com- 
ponent of the spectrum is constant. In this case, the superposition 
of this constant (in flux and shape) component and the contin- 
uum emission, of constant spectral slope and variable normali- 
sation, can result in spectral softening when the flux increases 
(see e.g. Taylor, Uttley, & McHardy 2003; Ponti et al. 2006; 
Miniutti et al. 2007). Furthermore, the presence of an absorber 
whose either the column density, ionisation state and/or cover- 
ing factor of the source varies, can also explain the observed 
spectral variations in the 2-10 keV band of the radio-quiet AGN 
(see e.g. Turner & Miller, 2009, for a review). Some cases of ex- 
treme absorption variations have indeed been reported so far in 
Seyf ert galaxies, with NGC 1365 being the most prominent case 
jRisaUti et al..2005; Risaliti et al..2007,: Risaliti et al..2009) . 

Clearly, the best option to avoid the effects of variable ab- 
sorption to the shape of the X-ray spectrum is to study the 
"hard" X-ray emission (i.e. X-rays at energies higher than ~ 20 
keV). At these energies, only neutral absorbers with a column 
higher than 10^"* cm"^ can affect the shape of the intrinsic spec- 
trum. Apart from these cases (which correspond to the so called 
"Compton" thick sources) hard X-ray light curves can in prin- 
ciple allow the study of the intrinsic continuum variability in 
AGN. AGN observations in hard X-rays have been performed 
the last twenty years by CGRO, BeppoSAX and INTEGRAL. The 
Swift/BAJi all-sky survey in particular has provided us with con- 
tinuously sampled light curves, which span periods as long as 
five years, for many AGN. 

Beckmann et al. (2007) have presented the results from a 
study of the first 9 month Swift/B AT light curves of 44 AGN. 
They found that ~ 30% of Seyferts exhibit significant hard X- 
ray variability on time scales of 20-150days, type 1 Seyferts are 
less variable than Seyferts 2, and a significant anti-correlation 
between luminosity and variability amplitude. More recently. 
Soldi et al. (2010) reported the results from a preliminary study 
of the flux variability of 36 AGN using data from the first 5 
years of ^wZ/f/BAT observations. Their results confirmed the 
hard X-ray variability - luminosity anti-correlation detected by 
Beckmann et al. (2007) at high energies. They also showed that 
an anti-correlation between variability amplitude and BH mass 
may also exist for Seyfert galaxies. 

In this work, we present the results from a variability study 
of the five brightest radio-quiet AGN in the recently published 
catalogue of Baumgartner et al. (2011), using the 5 years long 
BAT light curves that the same authors provide in the 20-50 and 
50-100 keV bands (the "soft" and "hard" bands, hereafter). Our 
main aim is to study the spectral variability of the sources with 
the use of "hardness ratios", i.e. by simply dividing the hard over 
the soft band light curves. Such ratios have been extensively used 
in the past for the study of the AGN spectral variability in the 2- 
10/20 keV band. Their biggest advantage is that they are entirely 
model-independent; if the hardness ratios are variable, then the 
spectral shape of the source has to be variable, irrespective of 
the underlying continuum spectrum, and of which spectral com- 
ponent is responsible of the observed variations. In addition, the 
presence (or absence) of a correlation between the hardness ra- 
tios and the source flux can indicate which model components 
vary (or not) in AGN. 

We restricted our study to the 5 brightest Seyferts in the cur- 
rent ^wZ/f/BAT catalogue because they are bright enough for an 
accurate estimation of their soft and hard band fluxes on time 
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Fig. 1. Swift/BAJ: light curves in the 50-100 keV (bottom panel), 
20-50 keV (middle panel), and full band (top panel) for the 
sources we studied. Solid circles and open squares indicate the 
2d (1 d in the case of NGC 4151), and the 20d-rebinned light 
curves, respectively. The time axis counts days since the start 
of the ^wZ/f/BAT observations. The dashed vertical lines in the 
NGC 4151 plot indicate the case of a flux variation whose evo- 
lution in the hard band is faster than the respective flux change 
in the soft band (see section 3.1 for details). 

scales as short as 20 days. As a result, we are able to use the 
hardness ratios to search for low-amplitude spectral variations 
on these time scales, almost continuously, over a period of 5 
years. This would not be possible to achieve with the study of 
energy spectra extracted over periods as short as ~ one month, 
due to low signal-to-noise ratio. In addition, the signal-to-noise 
ratio of even the 1-2 day binned soft band light curves of these 
objects is high enough to study their flux variations over a broad 
range of time scales, i.e. from years down to almost a day. 

The sample and the light curves we used are described in 
Section|2] Our results from the flux and spectral variability anal- 
ysis are reported in Sections[3]and|4l We discuss possible impli- 
cations of our results in Section |5] and we present our conclu- 
sions in Section|5] 

2. The sample and data analysis 

The Burst Alert Telescope (BA T; Barthelmy et al.2005) on-board 
the Swift (iGehrels et al.. 20041) is sensitive to X-ray photons in 
the 1 4-1 95 keV energy range. Baumgartner et al. (2011) released 
a catalog of sources detected in the first 58 months of BAT ob- 
servations. It consists of 1092 sources, detected at a significance 
level of at least 4.8cr. The majority of the sources in this cat- 
alogue are AGN (with 519 objects classified as Seyferts). We 
chose to study the 5 Seyferts which have the highest flux in 
this Swift BAT 58-Month Hard X-ray Survey catalogue, among 
all radio-quiet AGN. Source names and their hard X-ray fluxes 
are listed in Table [T] A summary of previous findings, regard- 
ing the hard X-ray emission of these objects, is presented in the 
Appendix. 

Baumgartner et al. (201 1) also provide light curves in eight 
energy bands: 14-20, 20-24, 24-35, 35-50, 50-75, 75-100, 
100-150, and 150-195 keV. These light curves are available 
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Table 1. The objects studied in this work. Average fluxes in the 14-195 keV band are taken from Baumgartner et al. (201 1). Columns 
3 and 4 list the BH mass measurement we use in this work, and the references for these values. The fractional variability amplitudes, 
as well as values when we fit the light curves with a constant (and number of degrees of freedom), are listed for the 20-50 keV 
band light curves (in parenthesis for the hard 50-100 keV light curves). They are calculated using light curves binned in 2 days (1 
day in the case of NGC 4151; see text for details). 



Name(Type) 


Flux 

(xlO""'ergs"'cm"^) 


F 

^ viir 


;r'/(d.o.f) 


BH„,,,, (Mo) 


Ref. 


NGC 4151(81.5) 


5.33 


0.370+0.005 


11058/1326 


4.6!^§xl0' 


1 






(0.37±0.01) 


(4272/1326) 




NGC 4945(82) 


3.01 


0.56±0.02 


2491/568 


1.4x10^ 


2 






(0.43+0.03) 


(1096/568) 






NGC 2110(82) 


2.97 


0.43+0.02 


2823/585 


2x10* 


3 






(0.45+0.03) 


(1256/585) 






IC 4329(81.2) 


2.90 


0.35+0.02 


1403/526 


1.3!i;OxlO'* 


4 






(0.32+0.05) 


(703/526) 




NGC 4388(82) 


2.76 


0.47+0.01 


2985/658 


6x10^ 


3 






(0.48+0.02) 


(1336/658) 







Notes. The references for the BH mass estimates are as follows: 1: Bentz et al. (2006), 2: Greenhill et al. (1997), 3: Woo & Urry (2002), 4: 
Mai-kowitz (2009) 
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Fig. 1. (Continued.) 



from HEASARCQ. They were extracted from the individual 
snapshot images from each ~5 min observation, and the reported 
count rates are corrected for off'-axis eff'ects. We added the count 
rate of the second, third and fourth band light curves to produce 
a combined light curve in the 20-50 ke V energy band (the "soft" 
band light curve hereafter; the error on the final count rate was 
calculated using the usual error propagation rules, i.e. Bevington 
1969). In a similar way, we produced 50-100 keV ("hard" band, 
hereafter), and the "full" band (i.e. 20-100 keV) light curves, 
by adding all the individual light curves in the respective energy 
ranges. 

The resulting light curves were then re-binned to 2 days and 
20 days (in the case of NGC 4151, we used a bin size of 1 day, 
as this is the brightest source in our sample). We produced these 
light curves by estimating the weighted mean (and its error) of 
all the points within each bin, provided there were at least 20 
points contributing to the estimation of the mean in it. The av- 



' http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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Fig. 1. (Continued.) 



erage "count rate over error" ratio (which is representative of 
the average "signal-to-noise ratio" of the light curve) of the hard 
band light curves is ~ 4 for NGC 4151 and -2-2.5 for the other 
objects. The ratio increases to ~ 7 and ~ 4 for the soft band light 
curves of NGC 4151 and the other objects, respectively. Given 
these values, the 2-d binned light curves at hand can reveal varia- 
tions which have an amplitude at least ~ 25% (50%) of the mean 
flux of the sources in the soft (hard) band. 

The same average ratio increases to 15 and ~ 5 - 7 for the 
20-d binned hard band light curve of NGC 4151 and of the other 
objects, respectively, and it is even higher (8 - 27) for the 20-d 
binned, soft band light curves. Given these values, the expected 
average signal-to-noise ratio for the hardness ratios is ~ 1 1 for 
NGC 4151 and larger than ~ 4.5 for all the other objects. These 
values imply that, even for the lowest signal-to-noise ratio light 
curves, we will be able to detect spectral variations which have 
an amplitude at least ~ 20% of the mean hardness ratio. 
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Fig. 1. (Continued.) 
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Fig. 1. (Continued.) 



3. Flux variability 

The soft, hard and total band light curves are plotted in Fig. [1] 
They cover ~5 years of almost continuous observations of the 
objects, and as we argued in the previous section, their length, 
sampling pattern and relatively high signal-to-noise ratio make 
them ideal for the study of the hard X-ray flux and spectral vari- 
ations of the sources. A visual inspection of the light curves indi- 
cates that all objects show significant variations. IC 4329 shows 
a gradual flux decay by a factor of ~ 2 in the first 2 years of 
the Swift monitoring (most pronounced in the total band light 
curve), and then its mean source flux level remains roughly con- 
stant with time. Other sources show larger amplitude variations 
on shorter time scales. For example, variations by a factor of 
~ 5, in less than year, can be easily seen in the light curves of 
NGC 4151 and NGC 4388. 

In order to quantify the presence of significant variations in 
the light curves we performed a simple test when fitting the 



data to a constant. The results are listed in the last column of 
Table 1, together with the number of the degrees of freedom 
(dof). The reduced values imply highly significant variations. 
The smallest value is 1.34 in the case of the hard band light 
curve of IC 4329, and, even in this case, the probability of the 
flux being constant is less than 4 x 10""*. To quantify the average 
variability amplitude of the sources we computed the fractional 
root mean square variability amplitude, Fvar, of the light curves. 
Following Vaughan et al (2003), this quantity is defined as fol- 
lows: 

where S ^ and x are the observed variance and mean of the light 
curve, and cf^j^ is the mean square error of all the measurement 
errors, (TeiT.i, of each light curve point, i.e. (ffgj^. - ('Zjf=i(^ln-i^/^- 

Since the mean square error, cr^„, is representative of the con- 
tribution of the measurement errors to the observed variance, 
S ^, Fvar should be an estimate of the intrinsic source variability 
amplitude over the time period sampled by the 5w/ff/BAT light 
curves, i.e. ~ 5 years. 

The results are listed in the 3rd column of Table 1 . On aver- 
age, the amplitude of the observed variations ranges from ~ 30% 
up to ~ 55% of the mean count rate. The average amplitude of 
all objects is comparable in both bands. This is not surprising, 
since the plots in FiglUshow that the same amplitude variations 
appear in the soft and hard band light curves of all objects. Fig.|2] 
shows a plot of the soft band Fvar values vs BH mass for the ob- 
jects in our sample (black, solid squares). The first thing to no- 
tice is that there appears an anti-correlation between variability 
amplitude and BH mass: smaller BH mass objects appear to be 
"more" variable. The linear correlation coefficient is -0.83, but 
due to the small number of objects, the probability of an intrin- 
sic correlation between Fvar and BH mass is just 0.08. The same 
anti-correlation between variability amplitude and BH mass has 
also been observed using both long term (e.g. Papadakis 2004) 
and short term X-ray light curves (e.g. O'Neill et al, 2005; Zhou 
et al., 2010) in the 2-10 keV band. 

The average variability amplitude of the AGN we studied 
is comparable to the average variability amplitude of the AGN 
long term, 2-10 keV light curves. For example, Sobolewska & 
Papadakis (2009) have recently studied the flux and variability 
properties of ten X-ray bright and radio-quiet AGN, using years 
long, 2-10 keV RXTE light curves. The open squares in Fig. |2] 
indicate their measurements. Clearly, at a given BH mass, the 
variability amplitude of the soft band Swift light curves is at least 
as large as that of the 2-10 keV light curves, despite the fact 
that the Sobolewska & Papadakis (2009) RXTE light curves were 
longer than ~ 7 years, i.e. longer than the Swift/Bat light curves 
we use in this work. 

3.1. Cross-correlation analysis 

We calculated the sample Cross Correlation Function, CCF(k), 
between the soft and hard band light curves following traditional 
methods (see e.g. Brinkmann et al., 2003). The CCF was com- 
puted at lags ^ = 0, + 1 , +2, +3 , . . ., days in the case of NGC 4151, 
while the lag bin size was equal to 2 days in the other objects. 
In all objects, the CCF peaks at zero lag with a maximum value 
which ranges between ~ 0.5 and ~ 0.7. At larger lag values the 
CCF decays to zero, with a similar rate in both the positive and 
negative lags. 

Given the excellent signal-to-noise of the NGC 4151 light 
curves, the resulting CCF for this object is much better defined 
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Fig. 2. 20-50 keV and 2-10 keV fractional variability ampli- 
tudes (black solid squares, this work, and open squares, from 
Sobolewska & Papadakis, 2009, respectively) plotted as a func- 
tion of BH mass. 
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Fig. 3. The NGC 4151 cross-correlation between soft and hard 
band light curves (open squares). The solid line (at lags < 0) 
indicates the CCF values at positive lags (see text for details). 



over a large range of time lags. The NGC 4151 CCF is plotted in 
Fig. E] (significant CCF peaks at positive lags in this plot means 
that the soft band variations are leading those in the hard band). 

The CCF peaks at a high value of ~ 0.8 at zero lag. The solid 
line at negative lags indicates the CCF{k > 0) values (i.e. it is a 
"mirror line" of the CCF function at positive lags). A comparison 
between this line and the CCF values at negative lags suggests an 
asymmetry towards negative lags in the sense that the CCF(fe < 
-10 days) values are larger than the CCF(A; > 10 days) values. 
If real, this would suggest that, on time scales longer than ~ 10 
days, the hard band lead the soft band variations. However, it is 
difficult to assess the significance of this result. 

We note though that there is at least one "clear" case where 
the hard band variations do appear to "lead" the variations in 



the soft band: the vertical dashed line in the middle and bottom 
panels of the NGC 4151 plot in Fig. [T] indicates the start of a 
flux drop "event" which appears simultaneously in both bands. 
However, although the hard band flux reaches its minimum flux 
level within ~ 40 - 60 days, the soft-band flux continues to de- 
crease for at least 60-80 days longer. This faster flux decline in 
the hard band is an example of a hard band "lead", which can 
explain the slight asymmetry towards negative lags we observe. 



3.2. Hardness ratio analysis 

We used the 20 d-rebinned soft and hard band BAT light curves 
(C7?2o-50keV and C/?5o_i()OkeV, respectively) to compute the hard- 
ness ratio, HR - C/?5o-iookev/C7? 20-50 keV, and study the spectral 
variability of the sources. We assume below that: a) a model fits 
the data well if the null hypothesis probability is larger than 1 %, 
and b) a model significantly improves the goodness of fit if the 
null hypothesis probability (i.e. "both models fit the data equally 
well") is less than 1%. 

First, we fitted the ffl?-time curves with a constant (here- 
after "model A"). This constant is basically equal to the mean 
HR value, HR. The resulting best-fit HR and;^'^ values, together 
with the number of degrees of freedom (dof), are listed in the 
second and third column of Table 2. These values indicate that 
the probability of a constant HR is less than 1% in all objects, 
except for NGC 2110 where the null hypothesis probability is 
10.1%. 

We then produced "colour-flux" diagrams, i.e. diagrams 
of the HR values plotted as a function of the total band flux 
(C7? 20-100 kev)- These diagrams are shown in Fig.|4](the full band 
count rates in these plots are normalized to the respective mean 
count rate). The NGC 4151 colour-flux diagram shows that the 
HR values cluster around ~ 0.4. The high value when we fit 
a constant to the NGC 4151 diagram is due to the presence of a 
few points, with small error bars, which scatter randomly around 
the mean HR. Although these points may indicate significant, 
intrinsic spectral variations which last for a short period of time 
(i.e. less than 20 days), the overall picture that emerges from the 
NGC 4151 diagram in Fig.|4]is that, despite the significant, large 
amplitude flux variations we observe in this object, the shape of 
the spectrum should remain roughly constant in this source. This 
is in agreement with the results form the recent work of Lubinski 
et al. (2010). 

On the other hand, the NGC 4388 colour-flux diagram (and 
perhaps the diagrams of NGC 4945 and IC 4329 as well) sug- 
gest flux related spectral variations: as the source flux increases, 
the HR values appear to decrease. This is identical to what is 
observed in the X-ray bright, radio-quite Seyfert galaxies in the 
2-10 keV band: the spectrum becomes "softer" with increasing 
source flux. In order to quantify the significance of this apparent 
trend, we fitted all colour-flux diagrams with a linear function of 
the form: HR = Q'+y6xC7?2o-iookeV (hereafter "model B"; during 
the fit we took into account the errors on both variables, follow- 
ing Press et al., 1994). The model B best fit results are listed in 
the last columns of Table 2, and the solid lines in Fig.|4]indicate 
the best-fit lines. 

In the case of NGC 2110, model B appears to improve the 
goodness of fit to the HR-fiux data of this source (Ax^ = 9 for 
1 dof; null hypothesis probability= 0.5%). However, as we men- 
tioned above, a constant HR model already provides an accept- 
able fit to the colour-flux diagram of this source. In addition, the 
model B best-fit slope is consistent with zero within the errors. 
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Table 2. Best-fit parameter values, values and number of degrees of freedom of Model A to the HR-time curves (2nd and 3rd 
column), and of Model B ° to the Hardness-Flux plots (see text for details). 



Name 


Model A 




Model B 








KR 


A''/(d.o.f) 




a 




NGC4151 


0.40 


215/97 


0.005±0.012 


0.391+0.003 


212/96 


NGC 4945 


0.52 


188/91 


-0.07+0.03 


0.461+0.011 


124/90 


NGC2110 


0.44 


115/97 


-0.014±0.023 


0.419+0.008 


106/96 


IC 4329 


0.31 


129/90 


-0.09±0.05 


0.303+0.007 


121/89 


NGC 4388 


0.43 


174/92 


-0.13+0.03 


0.418+0.008 


136/91 



Notes, (a) KR = a + fix (20-100 keV Flux) 



We therefore conclude that, just like NGC 4151, we do not find 
significant evidence for spectral slope variations in NGC 2110. 

Contrary to this, the best-fit slope value of NGC 4388 is dif- 
ferent than zero (the case of constant HR) at the 3.3cr level. In 
addition, according to the F-test model B provides a signifi- 
cantly better fit to the data (A^^ = 38 for 1 dof; null hypothesis 
probability^ 2.3x10"'*%) when compared to model A. We there- 
fore conclude that NGC 4388 shows significant spectral varia- 
tions, with a "softer when brighter" behaviour. 

The model B best-fit slope of NGC 4945 and IC 4329 is 
negative, and consistent with the best-fit slope of NGC 4388 
(within the errors). However, the negative slope in these objects 
is significant at just the 2.2cr and 1.8cr level. On the other hand, 
model B fits relatively well the colour-flux diagrams of these 
sources (null hypothesis probability^ 1.1% and 1.4%, respec- 
tively), while model A does not. In addition, the improvement 
of the model B best-fit, when compared to model A best fit, is 
highly significant in the case of NGC 4945 (A^^ - 64 for 1 dof; 
null hypothesis probability^ 1 x 10"^%). This is not the case in 
IC 4329. The colour-flux diagram of this source is broadly sim- 
ilar to the diagrams of NGC 4945 and NGC 4388, but there are 
quite a few points which cluster around a value of HR ~ 0.2, 
irrespective of their flux. We conclude that there are indications 
for flux related spectral variations in NGC 4945 (mainly) and 
IC 4329, but they are not as significant as in NGC 4388. 

In summary, NGC 4151 and NGC 21 10 do not show signifi- 
cant spectral variations, while we detect significant spectral vari- 
ability in NGC 4388. The situation is less clear in NGC 4945 and 
IC 4329: their colour-flux diagrams are similar to NGC 4388, 
but the presence of spectral variations is not highly significant. 

Since NGC 4388 is the lowest flux object in our sample, we 
investigated the possibility that the spectral variations we de- 
tect in this object may be due to unaccounted systematic un- 
certainties in its SwiftfBat light curves. For this purpose, we 
downloaded from HEASARC the Swift/BAT light curves of AX 
J1631. 9-4752, a source with a 20-100 keV flux similar to the 
NGC 4388 flux. This is an X-ray pulsar in a High Mass X-ray 
binary. Its spectrum is well fitted with an exponentially cut-off' 
power-law model, with a cut-off energy of ~ 10 keV. It remains 
constant (in shape) both during flaring and non-flaring periods 
(Rodriguez et al., 2006). Fig. |5] shows the colour-flux diagram 
for this source. A few HR values in this diagram are negative, 
because the source's flux is very low above 50 ke'V, due to the 
fact that its X-ray spectrum is very steep above above 20 ke'V. 
For the same reason, the HR values are on average smaller than 
the HR values of NGC 4388 and the other objects in our sample. 
The points along the x-axis indicate significant flux variations, 
whose amplitude is even larger than the amplitude of the X-ray 
variations we observe in NGC 4388. 



0.7 - 




(20-1 00 keV) Normalized Count R 



Fig. 5. Colour-flux diagram of AX J163 1 .9-4752. 

Despite these large amplitude variations, the colour-flux di- 
agram of this source appears to be remarkably flat. The solid line 
in the same figure indicate the model B best-fit line to the data. 
The best-fit slope is -0.0023+0.007, entirely consistent with the 
value of zero. We therefore conclude that the flux related spectral 
variations we detect in NGC 4388 are most probably intrinsic to 
the source, and are not caused by any systematic uncertainties. In 
fact, if the error on the model B best-fit slope is less than 0.01 for 
an object with no intrinsic spectral variations, and a count rate 
similar to (or less than) the count rate of the sources in our sam- 
ple (like AX J 163 1.9-4752 for example), then the model B best- 
fit slopes in the case of the NGC 4945 and IC 4329 colour-flux 
diagrams should be indicative of significant spectral variations 
in these sources as well. 

4. Spectral variability 

4.1. Interpretation oftlie "colour-flux" diagrams 

Although the hardness ratios can be used to detect spectral vari- 
ations in a model independent way, they cannot identify, un- 
ambiguously, which are the model components that are actu- 
ally responsible for the observed variations. The main reason is 
that in most cases (including our work) the hardness ratios are 
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Fig. 4. Hardness ratio (50-100 keV)/(20-50keV) values, estimated using the 20d-binned light curves, plotted as a function of the 
total (20-100) keV normalized count rate. The solid lines show the best-fitting linear function to the data. 



computed using light curves over energy bands which are quite 
broad, and as a result multiple components contribute to the ob- 
served count rate. In order to investigate in a quantitative way the 
constrains that the colour-flux diagrams in Fig. |4] can impose on 



current theoretical models, we performed the "experiment" we 
describe below. 

Thermal Comptonization is the most commonly accepted 
mechanism for the X-ray emission from radio-quiet AGN. We 
therefore considered a "Comptonization plus reflection" model. 
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similar to what has been used numerous times in the past to fit 
the high energy X-ray spectra of many AGN, and we computed 
the expected hardness ratios when the model parameters were al- 
lowed to vary. We produced "theoretical" colour-flux diagrams, 
which we then compared with the observed diagrams. 

We first used the th ermally comptonized continuum mode l 
Nthcomp in XSPEC (Zdziarski et al.. 19961: IZvcki et al..l999h . 
The main model parameters are the spectral slope, F, and the 
electron temperature, kT,,, which determines the high energy 
rollover of the spectrum. As for the other model parameters, we 
assumed a seed photon temperature of 10 eV, and a "diskblack- 
body" input spectrum (none of these two parameters affect sig- 
nificantly the shape of the spectrum at energies > 20 keV). 
We then added a reflection c omponent using the pexrav code 
dMagdziarz & Zdziarski.. 19951) . available in XSPEC, which cal- 
culates the reflected exponentially cutoff power law spectrum 
from neutral material. In all cases, we assumed an inclination 

angle of 45°, Fpexrav — rNthcomp, and -Ef.pexrav — ^7"e,Nthcomp- 

The reflection amplitude, R, was negative so that pexrav would 
output the reflection component only. We also considered neu- 
tral absorption with Nh - lO^-'cm"^ (similar to what has been 
measured in the past for the objects we studied - see Appendix), 
although it does not affect significantly the spectrum at the ener- 
gies we consider 

We downloaded the diagonal BAT survey instrument re- 
sponse from HEASARCE and we used XSPEC to simu- 
late SwiftfQP^ spectra for Y = 1.3,1.7,2.1. For each F, 
we considered three different electron temperatures: kT^ - 
50, 100, 200keV. For each pair of (Y,kTe) values, we used five 
different normalization values of Nthcomp, from A = 1 to 
A - 9 photon cm"^ s"', in steps of AA - 2 photon cm"^ 
s to produce 5 simulated spectra (our results do not depend 
on the particular choice of the model normalization, as long as 
the ratio (AmaxMmin) is large enough to reproduce the observed 
maximum-to-minimum flux variations). Finally, we added to 
the model a constant flux reflection componenfl The flux of the 
reflection component was such that /Jmean = 0.3, 1, and 1.5 in 
the case of the Nthcomp "mean" spectrum (i.e. the continuum 
when Amean - 5 photon cm"^ s"'). Obviously, in this scenario, R 
should decrease with increasing flux, since the flux of the reflec- 
tion component remains constant. 

For each model spectrum we computed the (20-100), (20- 
50) and (50-100) keV count rates, and the respective hardness 
ratio: (50-100 keV count rate)/(20-50 keV count rate). As a 
result, we were able to produce model colour-flux diagrams, 
which we plot in Fig. |6] The x-axis in these plots indicate the 
20-100 keV model count rate normalized to the same band count 
rate of the "mean" spectrum. The spread of the values along the 
X-axis is similar in the model and in the observed colour-flux 
diagrams (Fig.lUi. 

The results for the 7?,nean = 0.3, 1, and 1.5 case are plotted 
in the left, middle and right panel in this figure, respectively. 
Within each panel, the top, middle and bottom curves indicate 
the model colour-flux diagrams for F = 1.3, 1.7 and 2.1, respec- 
tively. Clearly, even the mean HR value of a source can provide 
information about its Rmean and F. For example, the expected HR 
value should be ~ 0.2 if /?niean =1.5 and F ~ 2.1, as opposed 
to ~ 0.55 - 0.7 in the case of Rtnemi - 0.3 and F ~ 1.3. This is 



^ http://heasarc.nasa.gov/docs/swift/results/bs58mon 
^ This choice corresponds to the case of reflection from a distant, 
neutral reflector, i.e. the putative molecular torus. If the torus is a few pc 
away from the central source, the resulting reflection spectrum should 
be constant on time scales of a few years. 



of course the case for the widest separated parameter values we 
considered, and the difference between the model colour-flux di- 
agrams becomes less pronounced for smaller differences in the 
input physical parameters. Nevertheless, a comparison between 
the observed and model colour-flux diagrams plotted in Figs.|4] 
and Fig.|6]can provide interesting results. 

For example, the fact that ///?ngc4151,ngc2iio ~ 0.4 implies 
that F < 1.7 and Rmew < 1 in both sources. If F were steeper 
than 1.7, then HR should be smaller than 0.4, while if /?n,ean ^ 1, 
then we should had observed a significant positive coiTelation 
between HR and flux (see the diagrams in Fig.|6]for F > 1 .7 and 
^mean ^1) which is not the case. 

Open squares, filled and open circles in each group of curves 
plotted in all panels of Fig. |6] correspond to kT^ - 50, 100 and 
200 keV, respectively. Clearly, the model HRs, do not depend sig- 
nificantly on kTe. At a given flux, HR decreases with decreasing 
kTe, as the 50 - 100 keV band flux decreases accordingly. These 
variations become more pronounced for F > 1 .7 and /?mean > 1 
(see the two bottom group of curves in the middle and right panel 
of Fig. 5). In this case, when kT^ changes from 200 keV to 50 
keV, HR decreases by ~ 0.1, which is comparable to the HR 
variation when F steepens from 1 .7 to 2. 1 (for the same kTe). 

In any case, the general trend of all the model diagrams in 
Fig. |6] is that of a positive correlation between HR and flux. 
Since the reflection component has a constant flux, and it is more 
pronounced at energies below 50 keV, the soft band count rate 
should remain roughly constant at low flux states; as the hard 
band flux will continue to decrease with decreasing flux, HR de- 
creases accordingly. And yet both the NGC 4 1 5 1 and NGC 2110, 
where HR remains constant, despite the significant flux variabil- 
ity of the sources. 

The filled squares and triangles in the NGC 4151 colour- 
flux panel (Fig. [Hi indicate the model colour-flux diagrams when 
F = 1.6 and 1.7, respectively, /?n,ean = 0.3, and kT^ decreases 
from 200 keV, in the lowest flux, to 50 keV for the highest flux 
spectrum. The agreement between the observed and the model 
colour-flux diagrams is quite good. In the case of NGC 4151, 
any model diagram with a constant kT^ will have a positive slope 
larger than the slope of the dashed and solid lines shown in this 
panel, and will obviously provide a worse "fit" to the data. In 
other words, a kT^, variation with flux is required in this source to 
reproduce its observed flat colour-flux diagram. In the presence 
of a constant flux reprocessed component, a power-law com- 
ponent which simply varies with flux does not result in a flat 
colour-flux diagram. The kT^, variation between 50 - 200 keV 
that our analysis implies is in excellent agreement with the re- 
sults reported by Lubiiiski et al. (2010) for the same source. A 
similar physical picture could also apply to with the observed 
colour-flux diagram of NGC 2110. 

A positive coiTelation is opposite to what we observe in 
NGC 4388 (and NGC 4945 and IC 4329 as well). Within 
the context of the model we considered, we can reproduce an 
///?-flux anti-correlation only if the intrinsic F correlates pos- 
itively with flux. Filled circles, connected with a solid line, in 
the observed colour-flux diagram of NGC 4388 in Fig. |4] indi- 
cate the model HR values when /?mean - 0.3, kT^ - 100 keV, 
and F increases from 1.3, in the lowest flux, to 2.1, for the high- 
est flux model spectrum. Although the model colour-flux dia- 
gram slope is steeper than the observed best-fit slope (a result 
which implies that the intrinsic AF in NGC 4388 is probably 
less than the AFmodei - 0.8 we considered), it is obvious that the 
observed spectral variations in NGC 4388 (and to a lesser ex- 
tent in NGC 4945 and IC 4329 as well) can be explained by an 
intrinsic spectral slope steepening with increasing flux. 
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Fig. 6. Plot of model colour-flux diagrams in the case of a "thermal Comptonisation plus a reflection component constant in flux" 
model. Different values for the (mean) reflection covering factor, power-law photon index and the temperature of the electrons in 
the corona have been considered: Rmean = 0.3, 1, 1.5, from left to right, F = 1.3, 1.7,2.1, from top to bottom (in each panel), and 
kTg = 50, 100, 200 keV, indicated with open squares, soUd, and open circles, respectively. 



5. Discussion and conclusions 

We used the 58-month long Swift/BAT monitoring light curves to 
study the flux and spectral variability of the five brightest Seyfert 
galaxies in the catalogue of Baumgartner et al. (201 1). Our main 
results are summarized below. 

- Both the 20-50 and 50-100 keV band light curves, of all five 
objects, are significantly variable on time scales as short as 
1-2 days. Their average variability amplitude is of the or- 
der of ~ 30 - 50% of the mean source flux. We also found 
evidence that the variability amplitude scales inversely pro- 
portional with the BH mass of the objects, a trend which has 
also been observed in AGN, but at energies below 10 keV. 

Since the column of the neutral material which obscures the 
X-ray source in these objects is less than 10^"^ cm"^, it cannot 
affect significantly the X-ray emission of the sources above 20 
keV. Consequently, the 5w(ff/BAT light curves can reveal the 
"true" variability behaviour of the central source in these, with- 
out the complications associated with the effects of the absorb- 
ing material at energies below ~ 10 keV. Therefore, the observed 
variability must be intrinsic to the central source. Furthermore, 
the similarity of the variability amplitudes at energies above 20 
keV and below 10 keV in AGN (for a given BH mass) implies 



that the main driver of the observed flux variations below ~ 20 
keV on time scales longer than 1-2 days may also be the intrinsic 
normalisation variations of the source continuum. 

- The soft (20-50 keV) and hard band (50-100 keV) hght 
curves, in all objects, are well correlated. We do not detect 
any delays down to 1 day, which is the smallest time scale 
we can sample with the light curves we used. In other words, 
similar variations appear, almost simultaneously, in both en- 
ergy bands. 

In the case of NGC 4151, which is the obj ect with the highest 
signal-to-noise ratio light curves, we detected a soft band delay, 
at time scales longer than ~ 10 days. This is the first time that 
such a delay within the X-ray band has been observed in an 
AGN, but it is not easy to quantify the significance of this result, 
and it is not clear what are its implications. If real, such a result 
implies that a physical mechanism (perhaps a change in the tem- 
perature of the hot plasma in the corona?) affects first the high 
energy part of the spectrum of the source, and then propagates 
to softer energies. 

- A hardness ratio analysis suggests that, despite the signif- 
icant flux variations, the shape of the 20-100 keV spec- 
trum remains constant in NGC 4151 and NGC 2110. On 
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the other hand, we detected significant spectral variations in 
NGC 4388: the hardness ratios decrease as the flux increases. 
A similar trend is also observed in NGC 4945 and IC 4329, 
but its significance is marginal. 

There is not a clear understanding of the physical parame- 
ters that drive the spectral evolution of AGN. As reported in the 
Appendix, previous studies of the sources in the sample have re- 
vealed hard X-ray spectral variations due to changes of either the 
photon index of the X-ray continuum, or of the reflection ampli- 
tude, and/or of the temperature of the Comptonizing electrons. 

The use of the (20-50) and (50-100) keV band light curves 
to compute the HRs means that we are not very sensitive in 
the detection of kTg variations, as long as kT^ > 50 keV. 
Interestingly though, the flatness of the colour-flux diagrams of 
NGC 4151 and NGC 2110 suggests that the electron tempera- 
ture should decrease with increasing flux in these objects. This 
in agreement with the results of Lubiiiski et al. (2010) in the case 
of NGC 4151. Furthermore, the average HR value in these ob- 
jects suggest that /?mean should be less than unity, and Fintrinsic 
less than ~ 1.7. 

Comparison of the colour-flux diagrams of NGC 4388 (and 
perhaps in NGC 4945 and IC 4329 as well) with the predictions 
of thermal Comptonization models supports the view that the 
main driver of spectral variability in this source is intrinsic spec- 
tral slope variations: Fintrinsic correlates positively with the source 
flux, while the reflection component flux remains constant. As a 
result, Robs should also decrease with increasing flux. 

It is difficult to identify the physical parameter responsible 
for the Fintrinsic-flux correlation that we observe in NGC 4388. 
An anti-correlation between kT,, and source flux, like the one 
observed in NGC 4151, could explain it, if the optical depth of 
the corona remains constant: in this case, a decrease in the tem- 
perature of the corona results in a steeper continuum. . On the 
other hand, the same NGC 4151 observations suggest that, as 
kTe varies, the optical depth of the corona also varies so that the 
Compton parameter y (and hence F as well) remains constant 
(for this source). Therefore, flux related optical depth variations 
could also explain the F- flux correlation, if kT^ remains con- 
stant. However, due to the fact that previous studies have indi- 
cated kTg variations in NGC 4388 (see Appendix), we believe 
that variations of the coronal temperature are the main driver of 
the spectral variability in this object. 

A "softer when brighter" behaviour is commonly detected 
in Seyferts at energies below 10 keV, both on short and long 
time scales. Sobolewska & Papadakis (2009) have found that 
the observed spectral variations of X-ray bright AGN in the 2- 
20 keV band can be explained if Finmnsic ^ ^2-iok v that 
were the case, the max-to-min flux ratio of ~ 7 in the case of 
NGC 4388 should imply a Fmax/Fmin ratio of ~ 1.2. The dashed 
line in the bottom panel of Fig.|4]indicates the expected colour- 
flux diagram when F = 1.3, for the lowest flux, and F = 1.7 
for the highest flux spectrum (so that Fmax/rmin ~ 1-3). The 
model colour-flux diagram has an amplitude which is higher 
than the amplitude of the observed diagram, most probably be- 
cause Rmean > 0.3 in this source. However, the slope of both 
the model and the observed colour-flux diagram are quite simi- 
lar. Therefore, our results, are in agreement and support those of 
Sobolewska & Papadakis (2009). 

The absence of significant spectral variations in NGC 4151 
and NGC 2110, although rare among Seyfers at lower energies, 
is not unique. For example, Sobolewska & Papadakis (2009) no- 
ticed that "NGC 5548 displayed limited spectral variations for 
its flux variability", and Papadakis, Reig & Nandra (2003) ob- 



served the same behaviour in PG 0804H-761. One important dif- 
ference between NGC 4151 and NGC 2110 and the other three 
sources in the sample is that they accrete at ~ 1% (or less) of 
the Eddington limit (see Appendix), as opposed to ~ 10% (or 
higher) for the other sources. Recently, Sobolewska et al. (201 1) 
studied the long term "spectral slope-flux" evolution of two well 
studied GBHs, namely GRO J1655-40 and GX 339-4. Their re- 
sults indicate that at accretion rates ~ 0.01 (or smaller) of the 
Eddington limit, the X-ray spectral slope is ~ 1.3 - 1.6, and it 
remains roughly constant over a large flux range (see the mid- 
dle panel of their Fig. 2; in particular the panel with the GRO 
J1655-40 data). This is exactly the case with NGC 4151 and 
NGC 2110 as well: they both have an accretion rate lower than 
~ 1%, a spectral slope flatter than ~ 1.7 (see e.g. Lubiiiski et al. 
2010, Winter et al. 2009) and do not exhibit significant spectral 
variations. These similarities raise the issue of different "spec- 
tral states" in AGN, just like in GBHs, with NGC 4151 and 
NGC 2110 being "hard-state" systems. 
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(iGreenhill et al..l9^). which rad iates at a rate of Lboi/Ledd = 
0.10 - 0.60 ("Mad eiski et al..200Cl) . Ginga observations revealed 
the presence of a hard X-ray source in this galaxy, which is heav- 
ily absorbed at energies below 10 keV by neutral material with 
an Nh as high as ~ 10^"^ cm"^ (Iwasawa et al., 1993). Large 
amplitude X-ray variations at energies above 10 keV, on time- 
scales less than <1 - 2 days (and as short as ~ 10"^ s) have 
been dete cted by RXTE, BeppoSAX and Su z aku observation s 
(Madeisk iet al..2000t iGuainazzi et al..2000l: lltoh et al..2008l) . 
These flux variations do not appear to be correlated with spec- 
tral variability as well (Guainazzi et al., 2000). Best-fit F val- 
ues range between 1.4 and 1.8 (see e.g. Iwasawa et al, 1993; 
Done et al, 1996; Madejski et al., 2000; Guainazzi et al., 2000; 
Beckmann et al., 2009; Winter et al., 2009). In addition, the high 
energy cut-off" in the X-ray spectrum is constrained to be in the 
range ~ 100 - 300 keV, while the reflection component appears 
to have a rather small value of R ~ 0.06 (Done et al, 2003). 

6.3. NGC2110 



6. APPENDIX 

We present below results from previous studies of the hard X- 
ray spectrum of the sources in the sample. 



6.1. NGC4151 

NGC 4151 is a nearby Seyfert 1.5 galaxy, and the bright- 
est persistent AGN in the 20-lOOkeV energy band (after the 
radio-loud blazar Cen A). It hosts a BH with a mass of 
4.6 "*^Q gXlO^ Mq, as es timated from reverberation mapping stud- 
ies (iBentz et al ,2006*), and radiates at a rate of of Lboi/Ledd - 
0.014 (Crenshaw & Kraemer,2007). 

Its spectrum has been extensively studied from radio wave- 
lengths to hard X-rays. Its energy spectrum above 20 keV is well 
fitted by a power-law model (PL, hereafter), with an exponential 
cut-ofF above energies 50-200 keV (see e.g. Piro et al., 1999 and 
2000; Beckmann et al., 2005 and 2009; Lubinski et al.,2010). 
Lubihski et al. (2010) in particular, have suggested that the cut- 
off" energy, E^, anti-correlates with the source flux, increasing 
from -5-80 keV, at high fluxes, to ~ 100 - 200 keV, when 
the source is at "dim flux state". The spectral slope of the PL 
component, T, is rather hard (compared to other. X-ray bright, 
Seyfert galaxies), with best-fit values being ~ 1.4 - 1.6 (see e.g. 
Zdziarski et al., 1996; Petrucci et al., 2001; Winter et al., 2009). 
de Rosa et al. (2007) have suggested that there may exist flux 
related, intrinsic spectral slope variations in this source. In addi- 
tion to the PL component, Ginga observations first established 
the presence of a "reflection hump" at energies above 10 keV 
(Zdziarski et al., 1996). This component has been detected by 
all major hard X-ray satellites since then, but with a variable 
amplitude, R, which ranges between 0.01 and almost 1 (i.e. the 
value expected in the case of a point source located on top of the 
accretion disc at relatively large height). Finally, it is well estab- 
lished that the X-ray spectrum of NGC 4151 is also affected by 
the presence of neutral material, perhaps with a variable cover- 
ing factor, and a column density, Nh, less than ~ 10^^ cm"^ (see 
e.g. Lubinski et al., 2010, and references therein). 



6.2. NGC 4945 

NGC 4945 is a nearby (3.7 Mpc; Mauersberger et al.l996) 
Seyfert 2 galaxy, with a BH mass estimate of a:1.4xlO''M0 



NGC 2110 is a nearby (z = 0.007579) Seyfert 2 galaxy. It 
hosts a black hole with a mass of 2x10^ Mq (Woo & Uitv,200J) 
and radiates at Lboi/Ledd - 1 0^^ - 10'^ ( | Evans et al.. 20071) . 
BeppoSAX observations (Malag uti et al.. 19991) revealed a steep 
power-law spectrum of F = 1 .9, and the presence of a neutral 
absorber with Nh~7x10^^ cm"^, partially covering the source. 
The same observations implied a rather low value for the re- 
flection covering factor (R<0.17 or < 0.5, if Ec = 1000 
or 50keV, respectivel y). The Second INTEGRAL catalogue 
(iBeckmann et al..2009b yielded the best-fit values of F = 2 and 
Nh = 4.3x10^^ cm"^ for the PL slope and the column den- 
sity of the absorbing neutral material, intri nsic to the source . 
Spectral studies of the Swift/BAT spectra (W inter et al.. 20091) 
provided the following estimates for the photon index and col- 
umn density of the neutral absorber: F - 1.54^^ ^^ and Nh = 
2.84+;!!^xl022cm-2. 

—0. 1 D 



6.4. IC 4329 

IC 4329 is a relatively nearby (z=0.0157; Wilson & Penston, 
1979), X-ray luminous AGN. It hosts a black hole with a mass of 
1.3+j,;gxlO'^Mo an d radiates at a ratio of Lboi/Ledd = 0.21+° °[; 
(lMarkowitz.20091) . BeppoSAX observations have shown signifi- 
cant variations in the 0. 1-100 keV flux emission from the source 
(Perola et al., 1999), which were not associated with signifi- 
cant spectral variations. COMPTON/OSSE observations have in- 
dicated the presence of a high energy cut-off in the hard X- 
ray spectrum of the s ource, at energies 250keV<E(;<1700keV 
jMadeiski et al..l995h . Subsequent INTEGRAL observations 
have indicated a smaller cut-off energy of ~ 80 keV (Beckmann 
et al., 2009). The best-fit X-ray spectral slope values at hard en- 
ergies range from ~ 1.4 to ~ 1.8 (see e.g. Beckman et al., 2009; 
Winter et al., 2009). There are indications of neutral absorption 
intrinsic to the source, with an Nh < 10^"* cm"^, and of the pres- 
ence of a reflection hump with a strength of R ~ 0.4 -1.2 (see 
e.g. Miyoshi et al., 1998; Done & Madejski, 2000). 



6.5. NGC 4388 

NGC 4388 is a nearby (z = 0.00842) nearly edge-on (;~78°) spi- 
ral galaxy hosting a Seyfer t 2 nucleus (Phillips & Malin, 1982 
iFilippenko & Sargent.l985b. It hosts a black hole with 



of 6x10'' M^ dWoo & Urrv.2002h and radiates at 
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LB„i/L Fnn ^ 0.1 (lElvis et al..2004 . SIGMA iPml etal..l99lh 
OSSE ('Johnson et al.JQQS") and BeppoSAX (iButler et al.. 19911: 
[Boella et al.,1991) observations during 13 years showed no 
spectral shape var iations of the P L. Subsequent INTEGRAL 
jBeckmann et al..2004.) and Suzaku dShirai et al..2008.) observa- 
tions revealed that the flux increased by a factor of 1.4 and 1.5, 
respectively (with a month to half-day timescales, respectively). 
The best-fit X-ray s pectral slope val ues at har d energies rang e 
from *1.3 to *1.8 dBeckmann et al.. 2Q09: W inter et al.. 20091) . 
with column density of Nh - 2 — 3x10^^ cm"^. The hard X- 
ray spectrum does not show a strong cut-off in the X-rays at 
£'<100keV, with the best estimate b eing Eg - 95^jTkeV, from 
INTEGRAL observations ( Heckmann et al.,2009). The reflection 
component was detected during the Suzaku observations and the 
reflection covering factors was constrained to be R = 1.3 - 1.5. 
There is barely a spectral shape variation between the high and 
the low-flux states, with the most important change being the 
normalization of the underlying power-law continuum from the 
nucleus. 
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